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Introduction. 

In vol 17, No. 3 (1905), of the ' Memoirs of the Imperial Academy of St. Petersburg/ 
M. Liapounoff has published an abstract of his work on figures of equilibrium of 
rotating liquid under the title " Sur un Probleme de Tchebychef." In this paper he 
explains how he has obtained a rigorous solution for the figure and stability of the 
pear-shaped figure, and he pronounces it to be unstable. In my paper in the 
' Philosophical Transactions 7# I had arrived at an opposite conclusion. 

The stability or instability depends, in fact, on the sign of a certain function which 
M. Liapounoff calls A, and which I denote A + 2 (B/) 2 /C/, where A is equal to 

a3[i(<r 8 ) 8 +2y-i(r 4 +sp,5]. 

M. Liapounoff tells us that, after having seen my conclusion he repeated all his 
computations and confirmed his former result. He attributes the disagreement 
between us to the fact that I have only computed portion of an infinite series, and 
only used approximate forms for the elliptic integrals in the several terms. He 
believes that the inclusion of the neglected residue of the infinite series would lead to 
an opposite conclusion. 

In my computation the function ^ 3 [|-(cr 2 ) 2 +2£ 4 ] — i^-or 4 is decisively negative, and 
being numerically greater than ^{(B/) 2 /C/ + [{, s]} 9 which is positive, the sum of the 
two is negative. The inclusion of the neglected residue undoubtedly tends to make 
this whole function positive, but after making the revision, explained in the present 
paper, it remains incredible, to me at least, that the neglected residue can amount to 
the total needed to invert the sign. 

It may be worth mentioning that in revising my work I notice that 
&3[sr(0"2) 3 +2£ 4 ] — -g-or 4 owes its negative sign to the term — ^cr 4 . This term arises 
from the energy of the double layer, called |DD. It comes from the portion of the 
term f 7rp 2 \e d ( 1 — Xe) dor , which gives rise to a term in e 4 with a negative sign. This 

* Series A, vol. 200, pp. 251-314. 
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term involves under the integral sign the factor |-(— + — J— 3G, all the other 

factors being positive. If we attribute to 6 and to <f> various values between ^rr and 
zero, we see that in part of the range the factor is positive and in other parts 
negative. A general inspection does not suffice to determine whether the positive 
portion outweighs the negative, as in fact it does. Therefore, in order to feel 
abundantly sure that no gross mistake had been made, I computed by quadratures 
the eight constituent integrals involved in the final result, and confirmed the correct- 
ness of the value found by the rigorous evaluation. 

The analysis of the investigation has been carefully examined throughout, and I 
have, besides, applied the same method to the investigation of Maclaurin's spheroid, 
where the solution can be verified by the known exact result. # 

As a further check, the formulae of the paper on the Pear-shaped Figure have been 
examined on the hypothesis that the ellipsoid of reference reduces to a sphere. The 
several terms correctly reproduce the analogous terms in the paper on Maclaurin's 
spheroid, but in effecting the comparison it is necessary to note that the variable r of 

1 — -g), whereas in the paper on Maclaubin's 

spheroid the corresponding variable r denotes 1( 1— -3), where r is radius vector and 

a the radius of the sphere. 

Dissent from so distinguished a mathematician as M. Liapounoff is not to be 
undertaken lightly, and I have, as explained, taken especial pains to ensure correct- 
ness. Having made my revision, and completed the computations as set forth here- 
after, I feel a conviction that the source of our disagreement will be found in some 
matter of principle, and not in the neglected residue of this series. I can now only 
express a hope that some one else will take up the question. 

In the revision of the computations, the methods now used are much better than 
the old ones. In as far as this paper is a mere repetition of the former work with 
improved methods, the results will only be stated in outline, but I now show how 
any of the ellipsoidal harmonic functions may be computed without approximation, 
and how the functions of the second kind may be found rigorously. 

The Cambridge University Press is now bringing out a collection of my mathe- 
matical papers, and when we come to the paper on the stability of the pear-shaped 
figure, the new methods of computation will be substituted for the old. 

This paper is supplementary to the former one on Stability, f and it will only be 
intelligible in connection therewith. As before, I refer to the papers in the 
' Philosophical Transactions 'J as "'Harmonics" and the " Pear-shaped Figure," 

* ( Amer. Math. Soc. Trans.,' 1903, vol. 4, p. 113, on "The Approximate Determination of the Form of 
Maclatjrin's Spheroid," and a further note on the same subject, recently sent to the same Society, 
t 'Phil. Trans./ A, vol. 200, p. 251. * 
I Yol. 197, A, p.- 461, and vol. 198, A, p. 301, 
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§ 1 . The Rigorous Expression for the Ellipsoidal Harmonic Functions. 

Rigorous forms have been found in previous papers for all the harmonics of orders 
up to the third inclusive. For harmonics of the fourth order rigorous algebraic forms 
may be obtained in all cases except when s = 0, 2, 4, but these are exactly the cases 
to be considered in this investigation. We have, then, to show how these functions of 
higher orders may be evaluated rigorously for an ellipsoid of known ellipticity. 

The only case required is that in which both i and s are even, and although all the 
forms might be evaluated by processes similar to those indicated below, I shall confine 
myself to this case. 

We have seen in " Harmonics" that if /3 denotes (1— k 2 )/(1 + k 2 ) of this paper, and ,u 
denotes sin 9, 

%' (fi) = P/ {fi) -/fy s+2 P/ +2 fa) + /% +4 P/ +4 ( J a) -...+ ( - f^WqiPi (p) 

-/3q s . 2 Pr 2 (fi) +P , q.- l P i - i (/*)"... + (~) |s /3 V; 0*). 
It is well known that 



and 



w-jr&v-w 



p . {fi) = (i+s)(i+s-l) (i-s + l) (1 _^ s [,,. (i-8)fy-8-l) ^_ s - 2 (1 _^ 



2 s . s ! 



l\2 2 (s+l) 

+ (*=£ ) (^'~ s ~ 1 ) (*'~*- 2 ) ( i-s-3) <-.-* ( -,_ 2) 
+ "~ 2!2 4 (s + l)(s+2) ~^ K ^ } 

Hence we may clearly write ^/ in the form 

f / (/*) = / sin 1 ' -f 2 sin 1 '" 2 cos 2 +f t sin' -4 cos 4 #-.... 

Since when s is not zero P/(l) = 0, and when s is zero P,(l) = 1, it follows that 
/o = (—f s fi is q . For the zonal harmonics (s = 0) this gives f a = 1. The determination 
of the other f's depends on that of the ^'s, which we shall consider later. 

Another form of ^/ (ju,) will be useful, viz. : 

$/ (fj) = a-b cos 2 0+c cos 4 0—d cos 6 0+e cos 8 0- .... 

It is obvious that 

a=/ 



5 



_ f i—2 f i(i—2) » 

c ~• /4+ 2m /3 + '2^^^ /o, 

rf -/« + 271]/* + — 2T2T— /2 + " 2 2 .3! /o> 
&c, &a, &c. 

Thus, when the f's are computed it is easy to obtain the a, b, c, e£, &a 

b 2 
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We know that (£/ (<j>) (the cosine function of <£) is the same function of — k' 2 sin 2 <j> 
that |3/ (/x) is of k 2 cos 2 0, except as regards a constant factor. 
Hence it follows that 



f K f2 K fi ft*' 6 

(£/(<£) = X a + 6_ g i n 2 ^ + c_ sm 4 <£H-(^sin 6 </>-f ... 

~. *^ "■ ** 

where X is a constant factor. 

Now I desire to define p/(ft) and C/(<p) exactly as in u Harmonics." 

This definition has already been adopted as regards fj/ (/&), but it remains to adjust 

the constant X so as to attain the same end as regards C/ (<£). 
When i and s are even, (£/ (<£) was defined thus : 



(£/ (<£) = cos s^> + j6p 4+2 cos (s + 2) <£ + /2 2 j3,+4 cos (5 + 4) <f> + . . . + /3 H * s) ^ cos i<f> 



4-/3p s _ 2 c °s (s—2)(f) + fi 2 p s ^ i cos (5 — 4) </> + ...+ /3 |s p . 



Since 



sm <£ 



2H. 



2r ! 



2 2r (r!) 2 2 ar " 1 (r-l)!(r+l)! 



cos2 ^+^=I 



2r I 



, f 



2 3r - 1 (r-2)!(r + 2)! 



cos 4<£— ... 



it follows that the term independent of <£ in (£,- (<£) is 



X 



2 



k 






1 



The term in cos 2</> in (£/ (<f>) is 



2X cos 2<£ 



JL J. 7 f( 

2*2 P 



' 3 , 1.3 2 

+ r .-C 



2.4 3 



K ' l,d*3 O t K 

jL . _ n „ _J_ 

r- 4 2.4.6 4 k 6 



The term in cos 4<]S> in <£/ (<£) is 



2X cos 4<j6 



l.o L , Zi K 

STI' S\4 ^ 



JL . o . -u . o 



2.4.6 4.5 



7 k /6 1.3.5.7 3.4 K fS 



w 6 



2.4.6.8 5.6 



K 



.8 



T* • • 



and so forth. 

In accordance with the definition to be adopted, these terms in the three cases 
respectively are : 1, cos 2</>, cos 4<£. Hence X must be chosen so as to fulfil that 
condition. 

Pursuing only the case of C t - ((/>) in detail, we have 



X 2 k 2 2.4 k 4 
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(£t (<£) = a'+b' sin 2 <j>+c' sin 4 cj>+ ... , 
we must have 

/2 ft 

a f = \a y V = X6 —=■ , c' = Xc —r , &e. 



2 ' 4 

K K 



Thus when / , f 2 , f±, &c, are found, it is easy to compute a, b 9 c, &e., and 
a', 6', c', &c. 

Our formulae tend to involve the differences between large numbers, and this defect 
becomes more pronounced as the order of harmonics increases. The fault is mitigated 
by using the forms 

p. fr) =/ sin* 6-f 2 sin*- 3 cos 2 0+ ... , 

Ct (<£) = a' + &' sin 2 <£ + c' sin 4 <£+.... 

In the case of a lower harmonic, however, such as the fourth, we may just as well 
use the form for ^ involving a, 6, c, &c, and powers of cos 2 6. 

We must now show how to complete the evaluation of the f's for the zonal 
harmonics. 

It appears, from p. 486 of " Harmonics," that, when i is even, we have to solve the 
equation 



ending with 



8a . = WiJiM i£iv8H*. 4 > 

H 4.1 2 +j8o— 4.2 2 +/8o— ... ' 

mkMh iz li , where {*,>} = (»+» (*-/+!). 

£ -j-pcr 



We are to take that root which vanishes when /} vanishes. 

Although the equation for fiar is of order i~l, yet at least for such an ellipsoid as 
we have to deal with, it is very easy to solve it by successive rapid approximations. 
It is clear that we may write the equation in the form 



(M+ 



1/Q2 

4 



^{♦,3^4} JF {i, 5} {i, 6} 
4.2 2 +/3o— 4.3 2 +y8o--.... 



fir = & {%, 1} [i, 2} . 



An analytical approximation is found by neglecting the continued fraction in the 
second term on the left, and we then obtain 

£cr = -2 + 2y/[l+$F(i-l)i(i+l)(i+2)]. 

If this value of j3cr is used in computing the first term of the continued fraction, 
and if the quadratic is solved again, we obtain a closer approximation. We then use 
the second approximation and include one more term in the continued fraction, and 
proceed until j3cr no longer changes. 



6 SIR G. H. DARWIN: FURTHER CONSIDERATION OF THE 

It is shown on pp. 486, 487 of " Harmonics" that 

g_2^ 1 jj8 a {^3}{^4} 

g 4.1 3 +j8cr— 4.2 2 +/3<r— ... 

2ff 4 ,_ 1 P 2 {j5}{i6} 

g 3 4.2 2 +/3o— 4.3 2 + / Bcr-... 

2(^ 1 ^{J7}{i8} 

q± 4.3 2 +i8cr— 4.4 3 + l 6<r— .. 



? 



It may be remarked that the factor 2 occurs in each of these equations on the left, 
excepting in the first one ; also we are to take q Q = 1. 

In the course of the successive approximations for the determination of /3cr, each 
of these fractions is naturally evaluated. Therefore it is only necessary to extract 
certain numerical values already found in the course of solving the equation for ficr. 

As a verification which shows whether the equation has been correctly solved, we 
have 

It is now obvious that we are able to find all the g ? s in terms of q 0i which is unity. 

1 — - /c 2 

We then multiply each q by its appropriate power of /3 or — — •$, that is to say, we 

1 + K 

form @ r q 2r for r = 1, 2, ..., ^i 9 and introduce the results into the formula for ^(ja). 

A closely analogous method enables us to find all the other types of function for an 
ellipsoid of known ellipticities, but, except for certain harmonics of the fourth order, 
it is not possible to obtain rigorous analytical solutions. Approximate analytical 
forms are given in " Harmonics," and the approximation may be carried further if 
desired. 

The following tables give the coefficients in the several functions for the critical 
Jacobian ellipsoid with which we are dealing :— 
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i. s. 


$. 


b. 


c. 


cl 


e. 


/ 


n o 

*2 2 


0-603374 
-0-039203 


0-923128 
0-923128 










4 

4 2 
4 4 


1 • 000000 
1*769147 
0-083965 


5-442161 
-36-154264 

~ 7-984389 


4-892138 
-44-93584 
95-562431 








6 
16 2 
t6 4 


1 • 000000 

- 8 • 4' 
3-78 


12-45814 
121-8 
-338-312 


29-55340 
439-425 
3680-303 


18-53561 
320-523 

4482-844 






8. 


1 • 00000 


23-29297 


103-90805 


155-9554 


74-7977 




10 


1 • 00000 


38-29978 


274-94458 


721-88640 


789-90216 


306-12784 




i s. 


a. 


//. 




d\ 


e. 


/'• 


n o 

*2 2 


0-603374 
- 0-039203 


0-076872 
0-076872 










4 
4 2 
4 4 


0-806905 
1-065020 
1-013640 


0-365661 
-1-812415 

-8-026680 


0-027371 

-0-187586 
8 • 000000 








6 
16 2 
16 4 


0-62544 
-1-1408 
1 • 0305 


0-64882 
1-5349 
-7-704 


0-12816 

0-4404 

6-944 


0-00669 

0-0264 

0-704 






8 


0-440664 


0-854891 


0-317488 


0-396793 


0-001585 




10 


0-289818 


0-924288 


0-552512 


0-120795 


0-011006 


0-000355 



* 



The second harmonics are here defined by fJ2 s (/*) = * 2 - q s 2 - k 2 cos 2 #, <ff 2 * = * 2 ~ ft 2 + k 2 sin 2 <j> (s = 0, 2) 
with q s 2 = | [1 + k 2 + a/(1 - /cV 2 )], with the upper sign for s = and the lower for s = 2. In the case of 
i = 4, 5 = 2, I had in the original paper inadvertently changed the sign both of fj 4 2 and $T 4 2 without, 
of course, introducing any error, since they occur always as a product, 
t These functions are only given in their approximate forms. 

As it is desirable to use the other form of 3f in the higher zonal harmonics, I give 
the coefficients f , f 2 , / 4 , &c, in these cases. It will be noticed how much smaller are 
the numbers involved. 



8 
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Coefficients of Terms in fj; (/*) when expressed in Sines and Cosines. 



i. s. 


u 


u 


u 


/e» 


U 


fio- 


6 

8 

10 


1*00000 
1 • 00000 
1 • 00000 


9-45814 

19-29797 
33-29978 


7-63713 

40-01416 

131-74546 


0-44035 

14-03323 

116-85134 


0*45235 
22-76398 


0-46728 



§ 2. The Rigorous Expression for the Harmonies of the Second Kind. 

The integral <&? denotes $/(z^) ©/(<>, and B/ denotes $f(v )-y- % s (v ). Thus 

3J/ is, in fact, the harmonic function of the second kind. B/ is clearly determinable 
from &\ 

In the original paper <U/ was found by quadrature, and this defect in my procedure 
is referred to by M. Liapounoff as a cause possibly contributory to the discrepancy 
between our results. Quadrature was not, perhaps, a very satisfactory method, and 
the defect will now be made good by finding these integrals in terms of the F and E 
elliptic integrals. It will appear that my former results were sufficiently near to the 
truth for practical purposes. 

The functions % s {v) or P/(^) are of eight types, determined by the oddness or 
evenness of i and s, and the association with a cosine or sine function of </>. In 
" Harmonics" the types are indicated by combinations in groups of three of the four 
letters E, O, O, S — denoting Even, Odd, Cosine, Sine ; for example, OES means i odd, 
s even, associated with a sine function. 

All the roots of the equation P/(V) or P/(p) = are real, and when the form of the 
function has been determined by the method of § 1 the equation may be solved. 
Hence these functions are expressible as the products of a number of factors ; and it 
is to be noted that it is not necessary to adopt the same definition as in " Harmonics," 
because the function may be multiplied by any constant factor, without affecting the 

For brevity, let 

n( K v~^) = ( ( cV- 2i 2 )(kV- ? /)...(kV- ? „ 2 ), 

1 

where k 2 is (l--/3)/(l+/3) of a Harmonics." The parameters, k and y, as elsewhere, 
define the form of the ellipsoid. 

An alternative notation will be needed, in which we write 



v 



ic sm 2 \j) 



, A/ = 1 — q/ sir i/f , A 2 = 1 — kt sir \ft. 
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9 



At the surface of the ellipsoid v = v^ \jj = y, and we shall, as before, write 
sin /3 = k sin y. At the surface of the ellipsoid we have then 



v = cosec )3, A x 2 = 1 — g/ sin 2 y , A 2 = cos 2 /3. 



In this notation 



w / A 2 



2 a 2 



tJkj ZA.2 • • i Za 

i'\sin a i/f/ sin 271 *// 



2 



A consideration of the eight types of harmonics shows that they may be written as 
follows : — 

Type. 



EEC, % n 2t (v) 



EES, P 2M+2 2J (.) 



ooc, p^r\v) 

OOS, $ 2n+1 2m (v) 
OEC, |k +1 » 



?i 



n 



II (kV-^ 2 ) 



cosec 2 " \jj II (A/), 



?i 



11 



= (kV-k 2 )* (kV- l) 1 IT (kV--^ 2 ) = cosec 2n+2 ^ A cos »// n (A 2 2 ), 



71 



ft 



(kV- 1)* n (kV-& 8 ) = cosec 2s+1 i|r cos ^ n (A/), 



7* 



» 



(kV-k 2 ) 1 n («V-g z 2 ) = cosec 2B+1 $ A n (A. T 2 ), 



?i 



Jt 



ki> II (K 2 v 2 —q x 2 ) 



cosec 2n+1 i/> n (A/), 



n 



n 



OES, P 2B+3 » =^(.V-K 2 )MKV-l)^n( K V-^) 



4 2« + 3 



EOC, P 2a+2 2{+1 (^) 
EOS, 1W» 



n 



kv (kV-i) 1 n (kV-^ 2 ) 



?i 



w (kV-k 2 ) 1 11 (#eV— <jr/) 



cosec" * " i// A cos i/i n (A/), 

1 



n 

cosec 2w+2 xjj cos t/f II (A/), 

i 



cosec 2w+2 i/r A n (A/). 

i 



Using p and © generically for any one of these and for the corresponding function 
of the other kind, we have 

» W € w = [» (v )T i o [f^7-(^-%(^-i/7p ' 

Or, changing the variable of integration to \jj, 



% = k [i (v )] 



c£i// 



■WW A 

To effect the integration the reciprocal of the square of ^ must be expressed in 
partial fractions. Inspection of the eight forms of functions shows that 



dp w: 



+ 



i ii ii 



JJ X 



f k 2 p 2 —k 2 g kV— 1 h kV ' 7 A x l(K 2 v 2 -q x 2 y 

with appropriate values of f 9 g, h, A X9 B x to be given hereafter. 
In every case but that of OES some or all of /, g 9 h are infinite. 



W-q* 2 j 
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C 
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In terms of t/» 






n 



WOT / a 



^ + ~ tan \jt+ T Biir.\b+ %-— 
g h i A x 



sin 4 xji B x sin 2 \fi 



A,' 



X 



hmk y 



This lias to be divided by A and integrated, and the result will be expressible in 
terms of the elliptic integrals 

J. o 



F{y) 



7 t 

aiji 



E{y) — A d\ji, n (y, A. c ) = 






Accordingly we require certain integrals, which are given on p. 313 of the "Pear- 
shaped Figure," but in somewhat different forms. Here and elsewhere k 2 denotes 
1 — k 2 , and (// denotes 1 — q x 2 . 

The integrals needed are as follows :— 



*y„:~2 



sm J xf; 







clxjj = 



1 in / \ , 1- -el / \ sin y cos y 

^ * (y) + -jj2 h (r) - - jtzz, a 



^ 



/c /c 



y tan 2 ^ ■ , . 

r - ' Cvtp 



1 



K 



/2 V / / r 



k 2 cos /3 
tan y cos /} 



k 



12 



> 



ry sin 2 i// 







A 



rfi/; = 



^ sin xb 7 1 



1 



o 



A,'A 



•a; 



2. 



1 1 

j * (y) - - 2 Hi (y) 

K K 

i n (y. a,) - i F (y) 



* * • » I i| 



y 



Sill ll/ 7 . 

— jUHi 

o A/A 



£ 



1 



L + _ 

2 ' .2 



2q x 2 \q f x 2 qx a 2 -?/ 



n (y> 4,) + 



1 



2& W* 5fV/ 



hF(y) 



1 



t^ / x sin y cob y cos ft 

o^ yl/ J "^ 2\ ^ XY) ~*~ 977~27T^ ^ 2\ A ™2 * 
^9* 9. x \ K —Qx ) Z< 1 x \K —^x ) &x 

The last two of these admit of considerable simplification, as will now be shown. 

It is proved in the " Pear-shaped Figure " that the II (y, A x ) elliptic integral 
disappears from the expression for !}fj(& for all cases up to the third harmonics 
inclusive. I have also proved numerically that for the Jacobian ellipsoid the like is 
true for all the even zonal harmonics up to the tenth inclusive, and for the tesseral 
harmonics i = 4, s = 2 and i = 4, s = 4 It is thus certainly true in all cases used 
by me, and I do not care to spend perhaps much time in proving algebraically the 
general truth of the law. 

The last two of our integrals will occur in the form 

f sin 4 ib j. D fv sin 2 xb 7 , 

** (1 ink •'J' JrouiwHi. *^ (/ mm* <** j4mB 

and I assume that the coefficient of II (y, A K ) in this expression always vanishes, as 

is proved to be true in all cases actually computed. 
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Hence 



1 1 1 
2B = •—-■ — — + ~— ~ 

x a f 2 a 2 K 2 -a 2 

1 x l x jL x 



Now the coefficient of F (y) in this same combination of integrals is 



1 



ZtQ ' x 



"1 1 

fix 9 a: 



In this we may substitute for J5 r its value, and thus find 



' y sin 4 xjj 

o A?A 



dty — B x 



ly sin 2 y 

o A r A 



d\jj 



1 pi / \ 1 -pi/ \ sin y cos y cos /3 

2<z,V-<z/) w 2 2 ,v,v-fc*) w 2?7R=^) A/ ' 



This expression together with (1) give all the required integrals, and it only 
remains to tabulate/, g, h, A x for the several types of function. 
For the sake of brevity I write 

C* = (?. 8 -ft?(?, 2 -ft 2 ) 2 ...(^ 2 -?» 2 ) 2 (x = 1, 2...n), 

the factor which would vanish being in each case omitted. 
When there is only one g, d is to be interpreted as unity. 

Table of Values of /, g, h, and A x . 



Type. 


* 

order of 
harmonic. 


,9 

rank of 
harmonic. 


/. 


(J> 


1 

- 

A. 

■ 00 


■**•& 


EEC 


2n- 


2^ 


00 


GO 


@x 


EES 


2n + 2 


2^ 


- k' 2 II (/<2 - q*f 

1 


n 

K'mq'J 
1 


GO 


42 (« 2 - <& 2 ) G x 


ooc 


2n + 1 


2tf+l 


GO 


1 


GO 


- i'2C x 


oos 


2n + l 


2^+1 


/c' 2 II (k 2 - ^2)2 
1 


GO 


GO 


- (* 2 - fe 2 ) a. 


OEO 


2n+l 


2t 


OO 


00 


n 

I% 4 
1 


( lx ^x 


OES 


2% + 3 


At- 


- k 2 k '2]I (k 2 - q*f 

i 


n 
1 


K 2 IT^4 
1 


?W (* 2 - qj) C x 


EOC 

! 

; 
1 


2% + 2 


2/+1 


GO 


1 


i 


_ A, 2,/ 2/7 
1/iB a x u x 


i 

; 

EOS 


2?i + 2 


2^+1 


1 


GO 


1 


- q2 (« 2 - q2) c x 



C 2 
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In all the types 



B r = 2 






1 



&-ti ' <i*- ( ii 



• » « | , 



1 



?/-?»'. 



with the omission of the term which would be infinite. 



We have generally 33 = 



a d 



- p (i>o). Hence by logarithmic differentiation of 



f (y ) dv c 

the expressions for the several types of p we find results given in the following 
"table i — -~ 

Table of the $5 Integrals. 



Type. 



EEC 



ET7IQ 



OOC 



OOS 



OEC 



OES 



EOC 



EOS 



order of harmonic. 



2w 



2w + 2 



2n + 1 



2n+l 



2n + 1. 



jQih "T" <0 



2n + 2 



2n + 2 



rank of harmonic. 



2t 



2t+l 



2t + l 






2/+ 1 



2^ + 1 



3&*-f 2 sin /?&;-'. 



1 



« 



|- sec 2 /? + 1 sec 2 y + 21/A^ 2 

i 



w 



| sec 2 y + 21/A^ 2 
i 



n 



£ see 2 /? + 21^2 

i 



?* 



| + 21/4,2 

1 



11 



■i + ^-sec 2 /3 + |sec 2 7 + 21/A^ 

i 



n 



| + 1 sec 2 y + 21 /A^ 



?i 



i + |~sec 2 ^ + 21/A a; 2 

i 



In the ease of the zonal harmonics (s = 0), (?*/* is always less than unity ; for 
harmonics of rank 2 one of the q x JK is greater than unity and the rest are less ; for 
rank 4 two of them are greater than unity and the rest less. 

For the zonal harmonics there is some gain in simplicity by putting sin 2 X = ^|. 
We then take the equation 

P&jO*) = a "~~b cos2 ^+ c cos4 #--... = o, 

and find all the n roots, say l9 2 ... n . 

If we solve the corresponding equation ^ 2n 2 (/x) = for the tesseral harmonic of 
rank 2, we find one root for cos 2 to be negative. If this root corresponds to 
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B u we must put cos* B x = 1 - £, so that q > = k» [1 + ( - cos* ft)]- Similarly for the 

harmonics of rank 4, two roots correspond with imaginary angles, and so forth. 
Subject to this explanation we may now regard the roots as defined by 6 U d 2 ...6 nt 
Since A/ = 1 — q 2 sin 2 y 9 we have A/ = 1 — ft 2 sin 2 y sin 2 6 X = 1 — sin 2 /3 sin 2 X9 and 

^(^) = cosec 2w yA/A/...A^ 
If we write 

D x = (sin 2 ^ - sin 2 .^) 2 (sin 2 # 7; - sin 2 # 2 ) 2 . . . (sin 2 9 X — sin 2 n ) 2 (n— 1 factors), 

our former C, may be written in the form k 4w ~ 4 D w and the several coefficients in the 
expression for gj/ may be expressed as trigonometrical functions— some of which may, 
however, be hyperbolic. 
We thus have 



2n* = ttt^i n (1 - sm2 fi sin 2 6 X ) 



2 



2k 



in- 



1 



w 



F(y)S 



1 



i jL/ r sin C7~ cos Uv 



x 



X 



n 



E ( y ) 2 






Dj, sin 2 6 X cos 2 #,, ( 1 — k 2 sin 2 # r ) 



?i 



+ k 2 sin y cos y cos /3 2 



1 



i D x cos 2 fr, (1 — k 2 sin 2 9 X ) ( 1 - sin 2 j3 sin 2 Z )_ 



This formula agrees with the result given for % 2 S (s = 0, 2) in § 4 of " The Pear- 
shaped Figure," although the formula is there expressed in terms of q 2 , and l/( K 2 — g 2 ) 
is replaced by its equivalent (1— 2<^ 2 )/(^y 2 . 

In the case of the even zonal harmonics of order i, all the ffs are real angles, and it 
facilitates the solution of the equation for 6 to note that, with rough approximation 
(improving as the order of harmonic increases), 



0i 



TV 



2t 



• ? 



e. 



3ir 
2?' 



{J t? . •—— ; „ , , , ? (7.1 t * 



3 



2i' 



~2» "' 



The following numerical values apply to the critical Jacobian ellipsoid : — 

For the fourth harmonic 6 X = 20° 15', 6 2 = 61° 11' ; the rough approximation gives 
22° 30' and 67° 30'. 

For the sixth 6 1 = 14° l'"9, 6 2 = 42° 12 /, 2, 6> 3 = 71° 8'*6 ; the rough approximation 
being 15°, 45°, 75°. 

For the eighth 6 X = 10° 43''1, # 2 = 32° 11''8, a = 53° 51''2, <9 4 = 76° 2l'"8 ; the 
approximation being 11° 45', 34° 15', 56° 45', 79° 15'. 

For the tenth Q x = 8° 40', B 2 = 26° 1', 3 = 43° 26', 4 = 61° 3', 5 = 79° 28'; the 
approximation being 9°, 27°, 45°, 63°, 81°. 



The values of the several %'s found by quadratures were in every ease too small ; 
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the correct values are given in the table below, I find that for ^ 4 quadratures gave 
too small a value by a -^jo^h part ; for f H 6 by a y-£otli P ar ^ '> ^ or ^ by an ~sVkh P^rt. 

The method which I have given above fails for the tenth zonal harmonic, unless 
we use logarithms of more than seven places ; and it is not worth while to 
undertake so heavy a piece of computation. I conclude by extrapolation that for gl 10 
quadrature (carried out on exactly the same plan as in all the other cases) gives too 
small a result by a y^-th part of itself I therefore augment in this case the result of 
the quadratures and find & 10 = 011640 ; this enables us also to compute 1S 10 . 

The following table gives the results of the whole computation" 76 : — 

Table of Logarithms of %? and $$/ Integrals. 



1, o. 


log E/ +10. 


log Ui». 


2 


9-6931231 


0929494 


2 2 


9- 


3330037 


4066504 


3 


9' 


54617 


20462 


4 


9' 


4332383 


2657402 


4 2 


9' 


24250 


39502 


4 4 


9' 


04753 


43121 


6 


9 


•2701270 


•3263106 


6 2 


9 


•14462 


39512 


6 4 


9 


•00632 


•42458 


8 


9 


•15835 


■35745 


10 


9-06595 


•36897 



§ 3. Note on § 15. " The Determination of certain Integrals: 

It has been found best to make some changes in this part of the work. 
The integrals to be evaluated are denoted 

*** sin 2 * 6 cos 2w e 



2p \ 2?i 
lx -2m 



r\ 2n 

11 2m 







Ai 2m A 



de, 






2 »of 



The integral 2p A 2m 2n may be made to depend on °A 2m 2n (which is the same as A 2m 
the original paper), and therefore I only evaluate the latter. 

These integrals were originally found as the differences of certain other functions, 
but it is not hard to give formulae for finding them directly. I have done this and 
recomputed the whole series of values, f 

* The only error of any moment which I have found in my previous work is that in some of the eases 
I had forgotten to introduce the factor #c in some of the gl's after effecting the quadratures. The error in 
my results from this oversight was fortunately not serious. 

f There is a misprint on p. 286. The function U$ should be 



QTT 



1 6 cos ft cos y 



[sec 4 /3 + sec 4 y + § sec 2 ft sec 2 y\ 
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The series of values °A 2n were also computed independently from the series given 
on p. 537 of " Harmonics/' 

A little consideration will show that the differences of the series of functions °A 2m 2n , 
with the signs of the odd differences changed, gives the series of functions 2p A 2m 2n . 
Stated analytically 

/ \ r \ r °A 2p—2r . 2r \ 2p—2r 

A table of the natural numbers °A 3m 2w is given, but the table of differences is not 
reproduced. 

The following is the result of the recomputation : — 

Table of the A Functions (Natural Numbers). 



n. 


°A ». 


°A 2 «. 


°A 4 '\ 


°A 6 ». 


2 


■7024906 8 


•034600 30 


53878 


132-38251 


2 


9505345 1 


4779482 2 


866414 


7*149917 


4 


•6559354 


8294118 1 


1590804 


1-8826818 


6 


5285432 


6160958 


7537022 


-9929042 


8 


•4543269 


•5084364 


5844941 


•6987687 


10 


4044492 


4419004 


4909693 


•5582507 


-^ a 


•3680175 


3958484 


•4306317 


•4755180 


14 \ 


•3399181 


•3616261 


•3878603 


•4203234 


16 


•3173978 


3349328 


3556021 




18 


•2988272 


•3133708 


•3295135 




20 


•2831734 








22 


•2697454 




- 





Table of Logarithms of the fl Functions. 



n. 


log a n . 


log i\\ 


log £l n 4. 


log tt 6 \ 





•2047610 1 


•0302912 


1-8667641 


2-7138144 


2 


9 


8993673 


•7715375 


1-6492558 


2-5319084 


4 


9 


•7729862 


6613000 


1-5528790 


2-4473549 


6 


9 


•6930884 


5897701 


1-4886134 


2-3893880 


8 


9 


6346685 


5365117 


1-4398970 


2-3446714 


10 


9 


5886267 


4939849 


1-4005014 


2-3080689 


12 


9 


5506357 


4585472 


1-3673618 


2-2769972 


14 


9 


5182995 


4281522 


1-3387296 


2-2499642 


16 


9 


4901531 


4015325 


1-3135070 




18 


9 


4652355 


3778480 


1-2909612- 




20 


9- 


4428427 


3565096 


1-2705764 




22 


9-4226147 









The results in the original paper were not so accurate as I had thought they were. 
There was a mistake in the differences which give the flt 6 series, affecting the values 
from n = 10 onwards, but as no use was made of the fl 8 series as published, the 
mistake did not affect the final result. 
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§ 4. Note on §§ 16, 17. The Integrals <r 2? cr 4 , £ 4 and &>/? />/? <£/• 

An improvement has been made in the method of computing all these. The 
functions to be integrated were written in every case with a common factor 
cosec 2 y(Ai 2 — r x 2 ) ; now this is equal to k 2 cos 2 6-Vk' 2 sin 2 <£. In consequence of the 
substitution of this value for the common factor we are able to obtain the result as 
the sum of, instead of the difference between, two numbers. Another consequence is 
that we can dispense with the series of functions denoted A_ 2 and fl_ 2 .' 

A single example of the way in which this change is applied will suffice. If we 
write 

/(A*, ) = a^-^A^+y %/-S% ^ 
/(A 2K 2 ) = a V-^A^+y^-S^ 



L 2w 
8 



>n = 1, 2, 



and denote by f(Q 2n ), /(0 2?l 2 ), corresponding functions with of, /3 ; , y\ 8 ; for a, /3, y, 8 
and O in place of °A ; it is easy to show that 

<r 2 = I cos 2 /3 cos 2 y {K 2 [/(A 2 2 )/(fl 4 °) + /(A 4 2 )/(0 2 »)-^(A 2 2 )/(0/)] 

+ k' 2 [/(V)./W) +/(A/)/(fl/)~ £/(A 2 °)/(0/)]}. 
The computations as revised gave 

cr 2 = '0136760, £ 4 = '000092343, cr 4 = '000176218. 

Taking log %$ = 9*5461687, I found 

^3[i(a- 2 ) 2 +2y~-lcr 4 - - '00050051. 

This differs by 4 in the seventh place of decimals from the old value. 
In evaluating co/, />/, <£/ when we use the form of $* involving f^f 2 ,f^ &c, we 
have to put 

[p 3 (/^)T = ^0 sin 6 0~F 2 sin 4 cos 2 0+ F 4 sin 2 cos 2 6, 

and jFq, -F 25 I?* are easily found from a, /3> y 3 8. 
I then write 

-£0 = /o^oj L 2 =s / F 2 +/2JP0J ^ = /0-P4 +/a^3 +f±Fo, &c. , 
and, when i denotes the order of the harmonic #/ concerned, write 

/(A 2n ) — £o f b A 2n °— i 3 * +4 A 2? , 2 + iy 4 l+2 A2 re 4 — ...(n = 1, 2), 
f(^2?ij = io t+6 A 2?i 2 -~Z'2 ? ' +4 A2 7i 4 'f -i 4 * +2 A 2w 6 — ,,,(n = 1 9 2). 
These functions are then combined to give the required integrals. 
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A similar notation enables us to evaluate <£/, which is given by 



H = 



so that 



6 



it 



6 



(k 2 cos 2 0+/c /2 sin 2 <£) &£ d6d<f>, 



# = - { K y(A 2 )/(n «)+K'y(v)/(o 2 )}, 

IT 



appropriate forms being attributed tof(A), f(Q>). 

The result of the ' procedure sketched is the following series of values, in which 
only the cases i = 6, ,9 = 2, 4, are derived from approximate forms. For the sake of 
comparison I add the approximate values of </>/ as computed from the formulae in 
"Harmonics." In the cases i = 2, s = 0, 2, the approximate results, derived from 
that paper, are multiplied by such factors as to make the approximate formulae for 
|J 2 (p) and (£ 2 (<£) agree with the exact one when [i = 1, <j> = 45° ; and for J3 2 2 (fi) and 
C 2 2 (<£) to make the coefficients of p 2 and cos 2 <£ agree with the exact formulae. 



Table of Logarithms of a)/, p-\ </>/. 















Approximate cj>f 


i 


6 # . 


log (Df +10. 


logpr'+lO. 


log </>f. 




from formula 
in " Harmonics." 










7-6310567 








2 





7 


•6714241 


7 


•0286816 


9-0051748- 


-10 


9 


•00518-10 


2 


2 


(-)5 


6818162 


(-)5 


0264000 


7-0397377- 


-10. 


7 


03981-10 


4 





8- 


0332932 


7 


3558076 


9*6886735- 


-10 


9 


68861 - 10 


4 


2 


(-)8 


25158 


(")7 


32157 


1-72739 




1 


72729 


4 


4 


8 


30779 


7- 


37092 


3-81610 




3- 


81612 


6 





7- 


96786 


7 


32449 


9-69177 - 


-10 


9 


69303-10 


6 


2 


(-)8- 


72778 


(-)*' 


94094 


— 




2- 


20562 


6 


4 


9- 


10094 


8 


13161 


— 




5- 


29999 


8 





7- 


78437 


6- 


96857 


9-75611 - 


-10 


9- 


76872-10 


10 





(-) 7-9838 


6-6024 


9-8473 


-10 


9-87800-10 



Note that w 10 is negative while </>io remains positive. 



The calculation of the integrals for i = 8 and i = 10 was very laborious, and as the 
results tend to present themselves as the differences betweeif large numbers, it is 
difficult to obtain accuracy with logarithms of only seven places of decimals. The 
integrals <f> are much the most troublesome ; indeed I do not claim close accuracy for 
<f> 8 ; and as it appeared to be impossible to compute <£ 10 to nearer than 10 per cent, 
from the formula, I computed the several constituent integrals for the tenth harmonic 
by quadratures and combined them to find <j> 10 . The results derived from the 
approximate formulae of " Harmonics " are given for the sake of comparison. They 
clearly give somewhat too large a value for the higher harmonics. I believe o> 10 and 
p 10 to be nearly correct. 

VOL. CCVIII. A. D 
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If allowance be made for the difference of definition adopted in this paper from that 
used in " Harmonics" as regards the second zonal harmonic, it will be found that 
ft) 2> o> 4 , co 6 , o) 8i o> 10 , when set out graphically, fall into an evenly flowing curve. The 
corresponding test for the p s is not quite so convincing, but there is nothing which 
implies a mistake. The values of p , p 2 , /> 4 , p Q fall well into line, and so do p 4 , p 6 , p 8 , p 10 , 
but there is a gentle elevation in the neighbourhood of p 6 . In consequence of this 
slight w^aviness of the curve I recomputed the ivhole again independently, after it had 
been recomputed and verified once, and special attention was paid to co and p 6 . 

§ 5. Final Synthesis of Numerical Results, and Conclusion, 

The several numerical values are combined just as in the original paper, but the 
numbers, of course, differ a little from those obtained before. The following table 
gives the final stage, inclusive of the additional terms now computed : — 



i. 


6*. 


(1) 
[i, s]. 


(2) 


(l) + (2). 


Bf/Cf. 


2 
2 
4 
4 
4 
6 
6 
6 
8 
10 



4 

2 
4 

2 
4 




+ •000138868 

•000000717 +■ 

•000092542 

•000001908 

•000000012 

•000031204 

•000003422 

•000000014 
+ •000012905 
-•000001030 ■•■ + 


•000219736 
000000970 
000154732 
000001190 

•ooooooooo 

'000031146 
•000002107 

•ooooooooo 

'000006671 
000007358 


- -000080868 
+ •000001687 
•000247274 
•000003098 
•000000012 
•000062350 
•000005529 
•000000014 
•000019576 
+ -000006328 


~~~ -12382 

- -08056 
+ •06273 

- -000355 
+ •0000017 
+ •019564 

- -000229 
+ -0000003 
+ -007505 

- -00667 


gls [1 (<r 


Sum = 
Numerator = 


•000265000 
- -000500513 


- -000235513 



I then find log D = 9*9840165, logZ= '6454565, log M= "9591963. From these 
we find c = L6 2 , 5 = jjf<fi 3 2 ; whence 



ra> 



-^2 -* 


- -0553908 


^2 2 


• 0008037 


6 = 

Denominator = 


•0316007 


- -0229864 
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The Numerator divided by the Denominator is — Sco 2 /47r/)e 2 , whence 



S- 2 
4:7rpe 

a> 2 



logJ^L_ 2 = (-) 8-01054. 



With = '141990, from § 7 of the " Pear-shaped Figure," we have 

27Tp 

co 2 +8a> 2 = w 2 [l-'1443066e 2 \ 

Thence we find 

f 2 = '195979c 2 , f£ = '603177c 2 . 

These values differ sensibly from the old ones. 

The moment of inertia with log a = 9*8559759 is given by 



^~^ r =J^[l + -157786c 2 ]. 
27rpk - 



The moment of momentum is 

3M 2 Z<o 

27rpk 



[1 + -085633c 2 ]. 



As before, we find the pear-shaped figure to be stable, because fche moment of 
momentum is greater than that of the critical Jacobian, provided that the infinite 
series does not amount to too great a sum. 

If e be the uncomputed residue of 2 « [>', s] + '-^/- 1, 1 then find, as before, that the 



moment of momentum is 

3Jf 2 aa> 



[l + *085633c 2 ~499*586ec 2 ]. 



2irpk 
The coefficient of e 2 will be positive and the pear stable, provided that 

499*5S6e < '085633, 

or 

e < -0001714. 

* 

The eighth zonal harmonic gave a contribution of '0000196, and the tenth of 
'0000063. These are respectively a ninth and a twenty-seventh of the critical total. 
The pear is then stable unless the residue of the apparently highly convergent series 
shall amount to more than 27 times the value of the last term computed. 
M. Liapounoff claims in effect to prove that this is the case, but to me it seems 
incredible. I look for the discrepancy between our conclusions in some other 
direction. 
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